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We calculate the Casimir interaction energy in d = 2 spatial dimensions between two (zero-width) 
mirrors, one flat, and the other slightly curved, upon which imperfect conductor boundary condi¬ 
tions are imposed for an Electromagnetic (EM) field. Our main result is a second-order Derivative 
Expansion (DE) approximation for the Casimir energy, which is studied in different interesting lim¬ 
its. In particular, we focus on the emergence of a non-analyticity beyond the leading-order term in 
the DE, when approaching the limit of perfectly-conducting mirrors. We also show that the system 
considered is equivalent to a dual one, consisting of a massless real scalar field satisfying imperfect 
Neumann conditions (on the very same boundaries). Therefore, the results obtained for the EM 
field hold also true for the scalar field model. 


I. INTRODUCTION 

The static Casimir force is a physical effect which man¬ 
ifests itself in systems consisting of a fluctuating (quan¬ 
tum, thermal,...) field in the presence of non trivial, 
time-independent boundary conditions [T|. The corre¬ 
sponding Casimir energy, may be characterized as a real¬ 
valued functional of function(s) which, under a certain 
parametrization, define the geometry of the boundaries. 
This observation can be used as the starting point for an 
approximation scheme, the Derivative Expansion (DE) 
originally proposed in Ref. [2] (for subsequent develop¬ 
ments see Ref. 0). The DE adopts its simplest form 
when the boundary conditions considered are ‘perfect’, 
i.e., they do not involve any parameter, and the geometry 
of the system is sufficiently simple, yet non trivial: One 
has two boundaries, one of them, R , is smoothly curved, 
and describable in terms of a single-valued ‘height func¬ 
tion’ ip, which measures the vertical distance of each one 
of its points to the flat boundary, L. In other words, 
such that R can be projected in terms of a single Monge 
patch, with L the projection plane. 

Under the assumptions above, one is clearly left with 
an energy which is a functional of a single function, ip, 
and also possibly a function of the parameters eventu¬ 
ally appearing in the definition of the boundary condi¬ 
tions (specially when they are imperfect). The DE is an 
approximation scheme for that functional, such that its 
leading-order term is tantamount to the proximity force 
approximation (PFA) [4]. 

The nature of the next-to-leading-order (NTLO) term, 
on the other hand, depends on the type of boundary con¬ 
dition being imposed on the field. For Dirichlet boundary 
conditions, it has been shown that it is always quadratic 
in the derivatives of the smooth functions ip, regardless 
of the number of spatial dimensions, d [5]. Therefore its 
contribution to the Casimir energy is the integral of a 
local function. The same happens for perfect Neumann 
conditions when d ^ 2, but the situation is qualitatively 


different when d = 2: the NTLO term becomes nonlo¬ 
cal in coordinate space, a phenomenon which, we have 
argued, is a manifestation of the existence of a massless 
excitation for the fluctuating field [5j. A similar effect 
may be seen to appear for the case of the EM field with 
perfect boundary conditions, also in d = 2. This is, as 
we shall show below, no coincidence, as both theories, 
scalar with Neumann conditions and EM field with per¬ 
fect boundary conditions are equivalent. 

In order to gain more insight about this issue, namely, 
the special nature of the NTLO contribution for the EM 
field in 2 + 1 spacetime dimensions with perfect bound¬ 
ary conditions, we perform here the following analysis: 
we consider a system with imperfect-conductor bound¬ 
ary conditions on both surfaces, and evaluate the leading 
and NTLO contributions to the DE. 

This study is of interest because of several reasons: on 
the one hand, one knows that the perfect-conductor con¬ 
dition is an approximation to a real, imperfect mirror. 
Besides, it will provide a way to cope with the infrared 
divergences which would appear for perfect conditions. 
Finally, note that, in spite of the fact that the system is 
defined in 2 + 1 dimensions, this analysis may be useful 
even when one considers the 3 + 1-dimensional case at fi¬ 
nite temperature. Indeed, thermal effects mean that one 
should take into account the contribution of the Matsub- 
ara modes |6j. Among them, the thermal zero mode be¬ 
haves as a d = 2 field and, as we have shown, an entirely 
analogous effect to the one has in d = 2 is induced [5]. 
The same can be said of a fluctuating Electromagnetic 
(EM) field with perfect boundary conditions at a finite 
temperature, since it can be shown to accommodate a 
Neumann like contribution. 

This paper is organized as follows: in Sect. ||l] we de¬ 
fine the system, starting by a description of the duality 
between the scalar and EM field models, and then con¬ 
structing the model for the EM field coupled to the two 
boundaries. We then define the respective effective ac¬ 
tion, functional of the shape of the deformed mirror. In 
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Sect. Ill we present results about the expansion of the 
effective action to second order in the departure with re¬ 
spect to the case of flat parallel mirrors. Then, in Sect. |IV| 
we deal with the DE for the Casimir energy, based on 
the results obtained in the previous Section. In Sect. [V] 
we present some examples where the NTLO correction is 
evaluated for imperfect Neumann boundary conditions 
in 2+1 dimensions. Sect. m contains our conclusions. 
Some technical details about the evaluation of the effec¬ 
tive action to the second order in the deformation are 
presented in an Appendix. 


II. THE SYSTEM AND ITS EFFECTIVE 
ACTION 

A. Scalar field / EM field duality 


hand, the Euclidean effective action which we shall calcu¬ 
late can be used (reinterpreted) as the high temperature 
limit of the free energy for a model in 3 + 1 dimensions, 
with boundaries defined by X 3 = 0 and X 3 = il>(x 1 , 2 : 2 ) 
(after a straightforward relabelling of the spacetime co¬ 
ordinates). 

The form of the boundary conditions imposed on the 
field at the boundaries (regarded as spacetime surfaces) 
is then 


d 2 <p(x) 


= 0 , d n p(x) = 0, (3) 

X 2 =0 X 2 =^(X«) 


where d n denotes the directional derivative along the di¬ 
rection defined by the unit normal to R, n^^y): 


S%-6£d a il)(: E||) 

^i+WKhW ' 


(4) 


Let us first see how a real scalar field in 2 + 1 dimen¬ 
sions with Neumann conditions may be described, alter¬ 
natively, in terms of an EM field with conductor bound¬ 
ary conditions. 

We first assume that we want to study a massless quan¬ 
tum real scalar field p(x), satisfying perfect Neumann 
boundary conditions on two static curves, denoted by L 
and R } the former assumed to be straight and the latter 
slightly curved. We use Euclidean conventions, such that 
x = ( 2 +), with fi = 0,1, 2, denote the (2 + l)-dimensional 
spacetime coordinates ( 2 + = imaginary time). Besides, 
we shall use the notation v» for a vector on a 1 + 1 di¬ 
mensional spacetime, e.g. un = (v a ) = (uojiq). We have 
introduced the convention, that we follow in the rest of 
this paper, that indices from the beginning of the Greek 
alphabet (a, /?,...) run over the values 0 and 1. No dis¬ 
tinction will me made between upper and lower indices, 
and their vertical position will only be decided having 
notational clarity in mind. 

The free Euclidean action for the vacuum field p is 
given by 


So = \ J d 3 x {dip ) 2 , (1) 

which is complemented by the assumption of Neumann 
boundary conditions on L and R. Regarding the Casimir 
energy calculation we just need static boundaries, but it 
is nevertheless useful to consider a more general, time- 
dependent expression for the curved boundary R. Thus, 
L and R are defined as the regions in spacetime satisfying 
the equations: 


L) x 2 = 0 , 

R) x 2 =if{x ll ) , (2) 

respectively. 

The reason to allow for such a time dependence is 
twofold; on the one hand the treatment of the problem 
is more symmetrical, and the physical case may still be 
recovered at the end by setting ij) = ip(x 1 ). On the other 


The scalar field p may then be mapped into the 3- 
potential Ap for an EM field, by means of the duality 
transformation: 


dpP <+ tiivp&i'Ap 


(5) 


where A M is a vector field. It is an immediate consistency 
condition of the above, by taking the divergence on both 
sides, that p is massless, □+ = 0, which we shall assume. 

Now, the boundary conditions © corresponds, via the 
duality transformation, to: 


^/5^ck^4/3 ( 2 ^11,0) 0 


[ nM (2j|) ^-pvp&vAp 


*2=0(*||) 


= 0 , 


( 6 ) 


for the Ap field. This may be expressed equivalently as 
the vanishing of the component of the EM field tensor 
which is ‘parallel’ to the respective surface; namely, the 
component of its dual (a pseudo-vector) parallel to the 
normal at each point vanishes. Let us consider that for R, 
since the situation for L may be obtained as a particular 
case, namely, if = 0: For R, introducing the projected 
component of the gauge field, Ao,( 2 :||): 


A*( ay) = e l f{x\\)Ap{x h 'il){x\\)) 

eaOll) = s a +S^d a if(x ll ) , (7) 


one gets on the surface the boundary condition: 


(2J11 ) 0 , (8) 

which for L simplifies to: 


£a@d a Ap(x\\, 0) =0, (9) 

i.e., the Ei component of the electric field vanishes for 
X 2 = 0. 

So, regarding the boundary conditions, we have a map¬ 
ping between Neumann and perfect conductor; for the re¬ 
spective free Euclidean Lagrangians, we note that, from 

©: 

-dppdpp = -^FpvFpv 


( 10 ) 
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so that the free scalar field action is mapped into the 
Maxwell action: 


\ J d 3 x(dip) 2 = SoM ^ S 0 [A] = \J d 3 xF 2 „(A) . 

(U) 

Now we want to deal with approximate boundary con¬ 
ditions, of such a kind that, for the scalar field would 
correspond to adding to the action an interaction term 
Sj localized on two mirrors, and containing the parame¬ 
ter /i, with the dimensions of a mass, such that perfect 
conditions are recovered when //—►(). More explicitly, 


SiM = ^ J d 3 x [ 6 (x 2 )(d 2 ip{x )) 2 


9(x\\)S(x 2 - ip(x\\))(d n (p(x))‘‘ 


( 12 ) 


where g(x y) = 1 + [dip(x \\)} 2 is the determinant of the 
induced metric on R, required to have reparametrization 
invariance. We use the same p, on both mirrors, since we 
will assume them to have identical properties, differing 
just in their position and geometry. 


B. The EM field model 

The approximate Neumann boundary conditions are 
then introduced in terms of the EM field, by adding to the 
EM field action the respective interaction term. Thus, we 
will work with the action: 


auxiliary fields are introduced in such a way that, if in¬ 
tegrated out, they reproduce the original action, S [A]. 
The corresponding equivalent action thus becomes: 

S(A,£ l ,£ r ) = 5 0 (A) +S 0 (£ l ,£r) - i f d 3 xJ l _ l (x)A fl (x) 

(18) 

with 


Jfj.(x) = S(x 2 ) 6 tia e a pdp€i,(x\\) 

+ S{x 2 - ip(x\\)) ^{x\\)e a pd^ K {x\\) , (19) 

and 


So(£l,£r) = ^ 


(£l 0||)) + v 9 { x \\) (£r(®||))" 


( 20 ) 


with g(x ||) = 1 + d a ipd a ip. 

The action <S(.A, £l,£r) ma y be seen to be invariant 
under gauge transformations, A^(x) —> A^{x) + d^ui(x) 
as a consequence of the fact that the ‘current’ is con¬ 
served. Since our next step amounts to integrating our 
the An, that action should be first given a gauge fixing. 
Gauge invariance assures the results are going to be in¬ 
dependent of the gauge-fixing adopted, thus our choice is 
dictated by simplicity. In this case that is the Feynman 
gauge, whereby one ads a gauge-fixing action S g f(A) to 
5o(A), to get the free gauge fixed action Sq(A): 


S 0 (A) S g (A) = S„(A) + IJ d 3 x(d-A ) 2 . (21) 


S(A) — <So(A) + iSi(A) (13) 

where 

S 0 (A) = \Jd 3 xF^F^ , (14) 

and 

5i(yl) = 5 l (A) + S r (A) . (15) 

where 

Sr04) = 1 1 ~ J d 2 x\i ^g{x\\)F a p{x\\)F a p{x\\) (16) 

with J" a ^(a;||) the EM field associated to „4 a (a:||), and /i 
the parameter introduced for the scalar field. The factor 
<Sl(^ 4) is defined by a similar expression, obtained by 
setting ip = 0: 

Sl(^4) = J d 2 x\\ F a p(x \\, 0 )F a p(x\\, 0) . (17) 

The interaction terms reproduce Eq. © when written in 
terms of the dual scalar field. 

Following standard procedures [7], we rewrite the ac¬ 
tion in an equivalent form, by using two auxiliary fields, 
and living on each one of the surfaces, to linearize 
the form of the terms localized on the mirrors. Those 


Now we define the effective action T[^>] by the func¬ 
tional integral: 

-Gdl _ JVAV£ l V£ r f d 3 x j^x)A„(x) 

6 ~ [VA e- 5 MA) ’ 

( 22 ) 

where the denominator has been introduced in order to 
get rid of one infinite factor which is irrelevant to our 
calculation: the effective action corresponding to the EM 
field in the vacuum, i.e., in the absence of mirrors. There 
are other factors we will get rid of in T, associated to 
the self-energies of the mirrors. These have the distinc¬ 
tive feature of being independent of the distance between 
the mirrors, and therefore they do not contribute to the 
Casimir force between them. 

We then integrate A, obtaining for T[^>] a formal ex¬ 
pression where we have to integrate over the auxiliary 
field, which are endowed with an action we denote by 
<Seff(£L,£R): 

e -r M = J X>£ l X>£ r (23) 

where 

Seff(£L,£R) = \ J d 3 xd 3 x'J^(x)G^ v (x,x')J fl (x') (24) 
= \J d 2 x\\d 3 x[ £ a (£||)A^b(£||,Z||)63(z||), 
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where 


then two-dimensional, the result being: 


G^(x,x') = 5^ 


d 3 k e ifc -(*-*') 

(27t) 3 k 2 


(25) 


is the A M propagator in the Feynman gauge, and we have 
introduced the four objects A"ab(®||,^J|)j where A and B 
adopt the values L and R. Their form is obtained by sub¬ 
stitution of the explicit form of in terms of the aux¬ 
iliary fields, performing integrations by parts, and using 
the respective 5 functions. Since we will use an expansion 
in powers of the deformation p, we only need them up to 
the second order in that expansion. See the Appendix 
for their explicit forms. 

We conclude this section by writing the effective action 
as follows: 


m = ^TrlogpK] , 


(26) 


where K is the 2 x 2 matrix of components defined by the 
kernels A'ab, while the trace operation acts over the A, 
B indices, as well as over the spacetime dependencies. 


III. EXPANSION OF r[ip] UP TO SECOND 
ORDER IN r? 


r (0) = tl- 
2 


d 2 k\\ 

(2tt) : 


log 


l-r 2 (|fc|||)e- 2 l fc " |o l . (30) 


with 


r(x) = 


x 

x + 2p ’ 


(31) 


while T and L denote the extent of the time and length 
dimensions of the system. We have extracted from 
an a-independent contribution. 

Thus, the energy density to this order has the form: 


£( 0 ) = 


47TC1 2 


dpp log 


l-( 


P 


p + 2 pa 


) 2 e~ 2p 


(32) 


This expression is well-defined for any value of pa; in 
particular, in the two limiting regimes corresponding to 
semitransparent mirrors, pa » 1 : 


£( 0 ) ~____ 

128tt p 2 a 4 ’ 

as well as for perfect mirrors, pa << 1 : 


(33) 


As we have already done in previous applications of the 
DE, we consider the effective action for an, in principle, 
time-dependent function i/j, taking the static limit at the 
end of the calculation. In that limit, the effective action 
becomes equal to the vacuum energy E, times T, the 
length of the time coordinates. 

Setting then tj)(x m) = a + p(x\\), we introduce the 
expansion of T in powers of p , up to the second order. 
Thus, 


,(o) „ _C(3) 1 

167 T a 2 


(34) 


This result corresponds, of course, to that of perfect Neu¬ 
mann boundary conditions. 


B. Next to leading order 


T(a, 77 ) = T ( 0 ) (a) + rW(a,77) + r< 2 >(a,jj) + ... (27) 


where the index denotes the order in p. The first order 
term can be made to vanish by a proper definition of a, 
and we consider the relevant, zeroth and second orders 
in the following subsections. They are: 

r (°) = iTr[logK (0) ] 

r (2) _ p(2,l) + p(2,2) ; (28) 


where: 

p(2.1) 


p(2,2) 


-Tr 

2 

-I Tl ' 


(K (0) ) 

(K (0) ) _1 K (1) (K (0) ) _1 K (1) 


(29) 


The second order term T^ 2 \ being quadratic in 77 , can 
be represented in Fourier space as: 

r<2> = \j (35) 

in terms of the kernel /( 2 h 

Collecting the two contributions to /^ presented in 
the Appendix, we obtain its full form. It may be rep¬ 
resented as the sum of two terms: fi , which is just 
quadratic in momentum and therefore gives rise to a lo¬ 
cal contribution to the effective action, plus another one, 
f n i where the dependence in k is inside the integrand 
of an integral over another momentum (pn), and gives a 
nonlocal contribution to the effective action; namely, 


A. Leading order 

The leading order term may be obtained rather 
straightforwardly, since the zeroth-order kernel is block- 
diagonal in momentum space, and the trace operation is 


with 


= k 


f {2) (k ||) = /,(fc„) + /„,(*„) ■ 

f d2 P\\ 


(2^(h,l)r 3 (b„l)^ 


(36) 


(37) 
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fnliH) = ~ J ^2^ lPll+k " la P(NI) 

+ r(|p|||)r(|p|| + fey |) + r 2 (| P || |)r 2 (|p|| + fcy |) 
- r 3 (|jj|||)r 2 (|jj|| + fcy|)e^ 2|Pi|+fe|||a | 


x S(|p|||)B(|p|| + fcy|) 


bu ■ (P\\ + fc ||)] 5 


where we have introduced the function: 

1 

X e 2ax — r 2 {x) 


(38) 


(39) 


It is quite straightforward to check that the previ¬ 
ous expressions render the proper limit for the perfect- 
conductor case, p —> 0 , under which r —> 1 : 

f i2) (k ||) /f (fc||) 


/n } (fc||) — 2 


rf 2 P|| _ 

(27r) 2 (e 2 l p nl a 


hi • (P|| + fc ||)] 2 

bulbil +fc|il 


l 

1)(1 - e - 2 lpn+*iik) 

(40) 


where also equals the kernel for a real scalar field 
with Neumann boundary conditions |7> . 


IV. DERIVATIVE EXPANSION 


The Casimir energy E, we have argued, is a functional 
of ip. Up to the second order, and recalling that ip de¬ 
pends on just one coordinate (xi), Ede, has the form: 


Ede [b] 


d.Ti 


V(b(xi)) + Z(b(a;i)) 


/ dip{x i) 
\ dx i 



(41) 

where V and Z are local functions of ip and /i. Taking 
into account the dimensions of the objects involved, we 
can write a more explicit form for V and Z: 


coOWO ^ = c 2 (pb) 

Mxi )} 2 ~ [bhi )] 2 ■ 


(42) 


where cq and c 2 , which determine the zeroth and second 
order terms, respectively, are dimensionless functions of 
their (also dimensionless) arguments. 

Regarding the cq coefficient, we find that: 


1 [°° 

co(fJ-a) = — dxx log 

47r J o 


1 - 


pa 


fe~ 2x 


■ (43) 


Note that, as shown in Fig. 1, co(pa) interpolates be¬ 
tween zero in the limit of large /id, and —C(3)/( 167 r) for 
pa —> 0, the result that corresponds to perfect Neumann 
boundary conditions (see Eq.(34|). 



FIG. 1. Coefficient 167rco/C(3) as a function of pa. 


On the other hand, the c 2 coefficient can be extracted 
from the kernel appearing in the second order term for 
the effective action, P^ 2 ), as follows: In the Taylor expan¬ 
sion around zero-momentum, we denote by a the coeffi¬ 
cient of term quadratic in the momentum: 

/ (2) (fc,|) = / (2) (0) + afcf + ... (44) 


Then, the c 2 coefficient is given by: c 2 (pa) = \ a . 

The contribution to c 2 coming from fi can be obtained 
from Eq. (37). Performing the angular integration, it 
reads 


pOO 

C 2 i{na) = 7 ^ / dx 
4tt Jo 


1 


(x + 2 pa) (x + 2 /ia) 2 e 2x — x 2 


(45) 

The calculation of the other contribution to c 2 , com¬ 
ing from f n i , is lengthy but straightforward. We expand 
the integrand in Eq. (38) in powers of fcy keeping just 
quadratic terms. Performing the angular integration we 
obtain 


C2 nil pa) = -— / dx -j 

4tt Jo (2pa + x) {x 2 — e 2x (2 pa + x) 2 ) 
xx 3 {e 2x x 4 (2pa + x) [8(pa) 2 (2a; 2 + x — l) 

+ 4pax(4x(x + 1) — 3) + x 2 (4x(2x + 3) — 5)] 

+ e 4x x 2 (2fj 1 a + x) [32(pa) 4 (x — 4)(2x — 1) 
+8(pa) 3 x(4x(4x — 15) + 23) 

+ 4(pa) 2 x 2 (28a; 2 - 94a; + 29) 

+ 2pax 3 (8x(3x — 8) + 19) + x 4 (4x(2x — 3) + 7)] 

— e 6x (2 pa + a;) 3 [—32(pa) 4 — 16(pa) 3 x — 12(pa) 2 x 2 
+ 12pax 3 + 3a; 4 ] + a; 7 } . (46) 


We have analyzed the behaviour of c 2 = c 2 j + c 2n ; 
both numerically and analytically. As expected, c 2 van¬ 
ishes as pa -A oo, which corresponds to the absence of 
mirrors. The most interesting limit is pa <C 1, since this 
case corresponds to almost-perfect Neumann boundary 
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FIG. 2. Coefficient C 2 as a function of [ia. 



that (up to the second order) DE approximation to the 
Casimir energy reads, for small ^a (i.e., close to the per¬ 
fect case): 


Eue [i>\ 


1 

167T 



[C(3) 


+ log[/*«*.)] (^) 


(48) 


and this constitutes one of our main results. The first 
term is the PFA for the Casimir energy (Epfa)- The sec¬ 
ond term contains the first nontrivial correction to PFA 
(F/ntlo) for an arbitrary boundary defined by ip. This 
equation shows that, although the DE is ill defined for 
Neumann boundary conditions en 2 + 1 dimensions, the 
non-analiticities in the DE appear only in the case of per¬ 
fect boundary conditions, that is, the parameter /i acts as 
an infrared regulator. The physical interpretation antic¬ 
ipated in Ref. [5] for the appearance of non-analiticities 
for Neumann boundary conditions is confirmed by this 
calculation. Indeed, for /r = 0 the field contain massless 
modes, that are not present for /i 7 ^ 0. Besides, note that 
the relation between the Casimir energies computed for 
different values of /i is encoded in the simple differential 
equation 


M- 


dE 


DE 


d/i 


16-7T 


dx 1 


ip{x 1) 2 L 


C(3) - 


under the assumption that: flip -C 1 . 


/<#(aq)\ i 

\ dx 1 y J 

(49) 


V. EXAMPLES 


FIG. 3. Logarithmic derivative — 1Q-k fia dc 2 ni /d^ia as a func¬ 
tion of fxa. 


conditions, a situation where the problems inherent to 
this case should start to manifest themselves. Indeed: to 
begin with, one can readily check that the integrand of 
Eq.( 46) goes like ~ \/x for /ra ~ 0, signalling the emer¬ 


gence of an infrared divergence, as anticipated in our pre¬ 
vious work. Moreover, it can be shown analytically that, 
when /ia <C 1 , 

c 2 (na) ~ c 2 nz(^a) - - 77 — log(A*o) . (47) 

l07T 


We have also checked this result by computing numeri¬ 
cally the logarithmic derivative of c 2 with respect to fxa, 
which indeed tends to — 1/(167 t) in the small-/ra limit. 
These results are illustrated in Figs. 2 and 3. Fig. 2 de¬ 
picts c 2 as a function of /ia, showing that it vanishes for 
large /id, while diverges in the opposite limit. As a quan¬ 
titative check for the small-^a behavior given in Eq.(47l, 


in Fig. 3 we show the plot of the logarithmic derivative 
/ia dc2ni/d^a. 

Collecting the results of this section, we can say 


Let us consider a function ip describing a parabolic 
boundary ip(x 1 ) = a + \x\/R facing a straight line, and 
approximate (n 7 ^ 0) boundary conditions (EM or Neu¬ 
mann, depending on the field considered). 

Note that a plays the role of the minimum distance be¬ 
tween the two boundaries, while R is the curvature radius 
of the parabola at its vertex. From Eq.([48|, we obtain 
the DE approximation to the Casimir energy, expected 
to be reliable, in this example, when e = a/R <C 1 . 

The zeroth-order (PFA) term, calculated from the first 
line in Eq. (48), reads 


A’pfa = — 


du 


C(3) 1 


C(3) f + °°_ = _ 

WnRJ ( e +i £) 2 Wy/lRei 


(50) 


The NTLO correction comes from the second term in 


Eq.(48). To evaluate the integral, we change the integra¬ 


tion variable to ip(xi) = za , so that: 

1 [2 f + °° \jz — 1 

Entlo = - sTr V 7 A * log {l,Ra} 


= ^I^V7 |1 + los(4t ‘ fle)1 - 


(51) 
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As expected, the NTLO correction to the DE expansion 
diverges logarithmically in the limit of perfect (/i = 0) 
Neumann boundary conditions, but is finite in the im¬ 
perfect case. 

Besides, we have found it noteworthy that, fixing the 
value of p and taking the limit e —> 0, the ratio between 
the NTLO and PFA terms becomes independent of p and 
reads 


-Entlo 2 

ge 


(52) 


This e log e behavior can be contrasted with the case of 
perfect Dirichlet boundary conditions; in Ref. J5] we have 
shown that, in 2 + 1 dimensions, 



4b M = - 


167T 


dx i 




[C(3) 


+ 


[1 + 6<(3)] / dip(x i) V 


12 


V dx\ J 


(53) 


FIG. 4. x ( |55| ) as a function of e = a/R (circle-line geometry) 
for pR = 0.1. x goes to 1 for small values of e, showing that 
the NTLO correction in this case is similar to the one of the 
parabola-line configuration. 


where the upper ( D ) denotes Dirichlet boundary condi¬ 
tions. Computing explicitly the integrals for the same 
example, we find 


4tlo _ (1 + ®C(3)) 
4*1 ” 24 C(3) ’ 


(54) 


which is linear in e, in contrast with the quasi-perfect 
Neumann case, that involves a logarithm. 

Let us now consider another example, that allows us to 
make contact with previous results in the literature [8]: 
a circle in front a straight line, which is the dimension- 
ally reduced version of the cylinder-plane case. Here, 
the function ip defining the curved contour is given by 
V’(aq) = a + R (1 — \J\ — x\/R 2 ), where R is the radius 


of the cylinder and a is, again, the minimum distance 
between the circle and the straight line. On physical 
grounds, we expect the results for either a circle or a 
parabola in front of a line to be very similar in the e<1 
limit. To check this assertion, since the integrals in Eq. 
(48) cannot be computed analytically for the circle, we 
have performed a numerical evaluation; in Fig|4]we plot 
the ratio: 


■Entlo 

* = (55) 

which compares the ratio between NTLO and PFA terms 
to the value it should have for the parabolic case, as a 
function of e. As expected, \ —> 1 for small values of 
e; exhibiting a similar behavior to the parabola-line case. 
Therefore, we conclude that the NTLO correction to PFA 
is proportional to e log e also for the circle-line geome¬ 
try, a result that was conjectured for perfect Neumann 
boundary conditions in Ref.(Hi. 


VI. CONCLUSIONS 

In this work, we present results which, we believe, may 
shed light on some of the properties of the DE approx¬ 
imation to the Casimir energy. In particular, we have 
studied a phenomenon pointed out in 5], namely, that 
for perfect Neumann boundary conditions in 2 + 1 di¬ 
mensions, the NTLO correction to the PFA cannot be 
written as the integral of a local term involving ip and 
up to two derivatives of that function. However, as the 
calculations presented here show explicitly, the NTLO is 
perfectly well-defined and local when the mirrors are im¬ 
perfect. In other words, the non-analyticity is an artifact 
of the idealization of the boundary conditions. And this 
is corrected by an imperfection, no matter how tiny. 

It is worth emphasizing that perhaps a taming of the 
non analyticity might also be obtained by introducing 
other, cruder infrared cutoffs, like a mass term for the 
field. However, as we have shown, it is sufficient to in¬ 
clude a rather mild and physically justified modification 
into the game, which consists of an imperfection in the 
Neumann conditions, parametrized by a constant that 
can be tuned to vary the mirrors’ properties. 

The same problem can be seen to appear when con¬ 
sidering the high temperature limit (for Neumann con¬ 
ditions) in 3 + 1 dimensions [5]. Mathematically, the 
integral in momentum space that defines Z{ip) has an 
infrared divergence. Based on analogies with results in 
the context of quantum field theory in non trivial back¬ 
grounds, we argued previously that the physical reason 
for the emergence of nonlocal corrections is the exis¬ 
tence of gapless modes, which happens only for Neu¬ 
mann boundary conditions. We have also argued there 
that, were that the case, for imperfect (and therefore 
more realistic) boundary conditions, the problem should 
be cured. We have shown here that this is indeed the 
case, by providing an explicit example. 














We considered the case example of the EM field in 
2 + 1 dimensions but, in what may be considered a by¬ 
product of our study, we have seen that it is dual to 
a real scalar field, in the understanding that perfect or 
imperfect conductor boundary conditions for the elec¬ 
tromagnetic field correspond, respectively, to perfect of 
imperfect Neumann boundary conditions for the scalar 
field. 

Regarding explicit results and examples, we have ob¬ 
tained the coefficients of the second order DE, depending 
on a parameter /r which measures the departure from per¬ 
fect boundary conditions, and applied them to evaluate 
the Casimir energy for the case of a parabola and a circle 
in front of a line. This enabled us to pinpoint the effect 
of the would-be infrared dominant contribution to the 
DE on the resulting energy, for specific geometries. We 
have also compared the results with those corresponding 
to Dirichlet boundary conditions, where the problem of 
non-analyticity of the NTLO correction is not present. 
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Appendix A: Intermediate results on the 
perturbative expansion for F 


We present here some technical details and interme¬ 
diate results corresponding to the calculation of the ef¬ 
fective action to the second order in the function 77 . We 
assume that ip(x\\) = a + r](xi 1 ), with a equalling the 
average of t/>. 

The term of order 1 vanishes, and the others can be 
written in terms of the expanded matrices K, of elements 
Kab (A, B = L, R), which have expressions that we 
present now. The zeroth order one 

K { 0 ) (z||,Z||) = K ( 0 ) (:r|| - :rf|) 

= J^e-ri^ mk] \ h (Al) 

where 

N+2/i/ 1 KN)e- |fc » |a \ 

2 VKN)e“ |fe " 10 1 ) ’ 

(A2) 

andr (+) = 

Regarding K^, we see that 

4r = 4l = 0 , (A3) 



and the off-diagonal elements are given by: 

^lr(*||>*||) = K RL( a! ||> a: ||) 

= _I / — jiH P ik ii'( x ii- x '\\) P -i fc iii° 

2 J (2tt) 2 

x [\ k \\\ 2 v{x\\) + ik a d a r](x ] \)\ . 


( 2 ) 

Finally, to the second order, = 0, and: 


+ 2 ) 

^RR 


( X \\, X '\\) = f 


d 2 h\ 




(2tt)2 

1 k a 

x + -\k^\ dad '^ 

- i\k\\Hx\\){d a v(x\\) +d' a r](x\ l )) 


(A4) 


(AS) 


On the other hand, both and are non-vanishing, 

but it may be seen that they do not contribute to the 
second order term under the assumption that the average 
of 77 vanishes. 

The Fourier transform of the inverse of (we need 
it to calculate the second order term) is given by: 


( 


K (0) ) '(fen) 

( 1 

V ~ r (\ k \\\) e 


2 


(|fc||| + 2 /r) (1 -r 2 
—r(|fc|| |)e - l fc n 1° 

|fc|||a ^ 


(N)e- 2|fc » |a ) 


(A 6 ) 


The second order term, T^ receives two contr ibu tions, 
which we have denoted by T^ 2,1 ^ and T^ 2,2 ^ in (28 1 . For 
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each one of them we introduce a momentum kernel, and 

r(2,a) = lj |^/ (2 ’ a) (fc ||)|r?(fc,|)| 1 2 , a = 1,2. (A7) 


An explicit evaluation of those two kernels yields: 


/ (2 ’ 1} (fc ||) 



d 2 P || r 2 (| P |||) 

(2tt) 2 e 2 “IP||l_ r 2(| p|| |) 


p 

IpiiI+ 2 p 


r ^ 2 P|| r 3 (|p|||) [P||-(P||+fc||)] 2 

J (2tt) 2 e 2 “IP||l_ r 2(| p|| |) |P||||P||+fc||| ’ 


(A8) 


,(2,2 ) (k , = _ [ d2 P\\ j_ r 2 (| P |||)r 2 (|p|| +fe|||) _ 

" J (27r) 2 | (e 2a l p nl — r 2 (|p|| I)) (e 2a l p ii +fc »l — r 2 (|p|| + fcy I)) 

_ KbiilMbu + fe|il) _) bn • On + fc n)] 2 

+ (e 2 “l p »l-r 2 (|p|||))(l^r 2 (b||+fc|||)e- 2a l p »+ fc nl) J |P||||P|| + fc||| 


(A9) 
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